Remarks 

This is in response to the final Office Action mailed on October 19, 2006. Claims 1 and 
23 are amended, support for the amendments being found, for example, at page 7, lines 23-26 
and Figures 2 and 4 of the present application. Claims 1-9, 12-31, and 34-45 remain pending. 
Reconsideration and allowance are respectfully requested for at least the following reasons. 

I. Specification 

The specification is objected to as failing to provide proper antecedent basis for the 
claimed subject matter. The correctness of the objection is not conceded. However, the claims 
are amended to address this issue. Specifically, the claims are amended to remove the recitation 
of a "computer-readable medium." Reconsideration and removal of the objection to the 
specification are therefore requested. 

II. Claim Rejections - 35 U.S.C. 101 and 112 

In sections 5-11 of the Office Action, the Examiner rejected claims 23-31 and 34-45 
under 35 U.S.C. §§ 101 and 112. The correctness of these rejections are not conceded. 

With respect to sections 5-9, claim 23 has been amended to remove reference to a 
"computer-readable medium." Claim 23 is now directed to a tangible optical watermark for use 
on a document. It is respectfully submitted that amended claim 23, and for the same reasons 
dependent claims 24-31 and 34-45, comply with 35 U.S.C. §§ 101 and 112. 

Additionally, claims 1 and 23 have been amended in response to sections 10 and 11 the 
Action, whereat the Examiner objected to claims 1-9, 12-31, and 34-45 under section 112. 
Specifically, claims 1 and 23 have been amended to make clear that the decoder matches the 
phase unmodulated (original) dot pattern. 

Reconsideration and allowance of claims 1-9, 12-31, and 34-45 are therefore requested. 

III. Claim Rejections - 35 U.S.C. § 103 

In sections 12-19 of the Action, the Examiner rejected claims 1-4, 6-9, 13-16, 20-26, 28- 
31, 35-38 and 43-45; 5 and 27; 12 and 34; 17 and 39; 18 and 40; 19; and 41 and 42, respectively, 
under 35 U.S.C. § 103(a) over the combination of Koltai et al., U.S. Patent No. 6,104,812 
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(Koltai) in view of Wang, U.S. Patent No. 6,252,971 (Wang) or the combination of Koltai and 
Wang in view of various tertiary references. These rejections are respectfully traversed, and the 
correctness of these rejections is not conceded. Reconsideration is requested for at least the 
following reasons. 

In section 12, the Examiner rejected claims 1-4, 6-9, 13-16, 20-26, 28-31, 35-38 and 43- 
45 under 35 U.S.C. § 103(a) as being unpatentable over Koltai in view of Wang. As 
acknowledged by the Examiner, Koltai fails to disclose or suggest embedding each latent image 
object into its respective watermark layer by phase modulation. The Examiner stated that Wang 
teaches using phase-modulated stoclustic screens (column 2, lines 32-35), and that one of 
ordinary skill in the art would have been motivated to modify Koltai to specifically embed each 
latent image object into its respective watermark layer by phase modulation. 

The phase modulation taught in Wang is based on stoclustic screens. The stoclustic 
screens disclosed by Wang do not have a fixed frequency and angle. The decoding process in 
Wang is therefore by stoclustic screens, which is performed by correlating un-modulated cluster 
positions with the modulated cluster positions. 

There is no teaching in Wang as to how one skilled in the art could simply apply 
stoclustic screens in watermark layers of fixed frequency (AM screen), as suggested in the 
Action. Furthermore, there is no motivation to apply stoclustic screens in multiple superposed 
watermark layers, since the process would have been extremely difficult and not obvious for one 
skilled in the art to simply apply stoclustic screens in watermark applications. 

In contrast, claim 1 has been amended to more clearly define over the combination of 
Koltai and Wang. Specifically, claim 1 has been amended to include the limitation that the dot 
pattern has a fixed frequency for each direction in the dot pattern. This configuration correlates 
with the embodiments discussed in the description and drawings, for example, on page 7 and 
Figures 2 and 4, corresponding to amplitude modulated (AM) screens. In the invention as 
defined by claim 1, decoding of the optical watermark requires demodulating the modulated dot 
pattern with a phase-modulated, fixed frequency dot pattern. 

It is respectfully submitted that Wang does not disclose or teach a dot pattern with such a 
fixed frequency. Wang does mention decoding with a "checkerboard pattern" (column 4, lines 
36-44 and FIG. 7); however, the purpose of the checkerboard pattern of Figure 7 is to simplify, 
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for illustrative purposes, the description of stochastic clusters. Actual stochastic clusters are not 
positioned like a checkerboard, as evidenced and shown in Figure 6 of the attached paper 
"Stochastic Clustered-Dot Dithering" by V. Ostromoukhov et al., Journal of Electronic Imaging , 
October 1999, Volume 8, Issue 4, pp. 439-445, published by SPIE. See the Appendix. 

It is therefore respectfully suggested that Koltai in combination with Wang fails to 
disclose or suggest embedding each latent image object into its respective watermark layer by 
phase modulation, where the dot pattern has a fixed frequency for each direction in the dot 
pattern. Reconsideration and allowance of claim 1, as well as claims 2-4, 6-9, 13-16, 20-22 that 
depend therefrom, are respectfully requested for at least these reasons. 

Claim 23 is directed to a tangible optical watermark for use on a document. Claim 23 
recites, among other limitations, at least one latent image object embedded into each watermark 
layer by phase modulation, where the dot pattern has a fixed frequency for each direction in the 
dot pattern. Claim 23 is allowable for at least reasons similar to those provided above with 
respect to claim 1. Reconsideration and allowance of claim 23, as well as claims 24-26, 28-31, 
35-38, and 43-45 that depend therefrom, are therefore respectfully requested. 

Regarding sections 14 -19 of the Office Action, claims 5 and 27; 12 and 34; 17 and 39; 
18 and 40; 19; and 41 and 42, all depend respectively from one of claims 1 and 23. None of the 
secondary references remedies the shortcomings of Koltai and Wang noted above. Since the 
Examiner's objections in sections 14-19 are made in view of the combination of Koltai and 
Wang in view of various tertiary references, it is respectfully submitted that claims 5, 12, 17-19, 
27, 34, and 39-42 are allowable for at least the same reasons as those provided above with 
respect to claims 1 and 23. Reconsideration and allowance are therefore respectfully requested. 
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IV. Conclusion 

Favorable reconsideration in the form of a Notice of Allowance is respectfully requested. 
Please contact the undersigned attorney with any questions regarding this application. 



Respectfully submitted, 
MERCHANT & GOULD P.C. 
P.O. Box 2903 

Minneapolis, Minnesota 55402-0903 
(612) 332-5300 



Date: March 19, 2007 /Robert A. Kalinsky/ 

Name: Robert A. Kalinsky 
Reg. No.: 50,471 

Attachment: Appendix 
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APPENDIX 

Attached hereto is the document "Stochastic Clustered-Dot Dithering" by V. 
Ostromoukhov et al., Journal of Electronic Imaging , October 1999, Volume 8, Issue 4, pp. 439- 
445, published by SPIE. 
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STOCHASTIC CLUSTERED-DOT DITHERING 



Victor Ostromoukhov 1 , Roger D. Hersch 2 
Ecole Polytechnique Federale, Lausanne (EPFL) 

Abstract 

A new technique for building stochastic clustered-dot screens is being proposed. A large dither matrix comprising 
thousands of stochastically laid out screen dots is constructed by first laying out the screen dot centers. Screen dot 
centers are obtained by placing discrete disks of a chosen radius at free cell locations when traversing the dither 
array cells according to either a discretely rotated Hilbert space-filling curve or a random space-filling curve. After 
Delauney triangulation of the screen dot centers, the maximal surface of each screen dot is computed and iso-inten- 
sity regions are created. This iso-intensity map is converted into an anti-aliased grayscale image, i.e. into an array 
of preliminary threshold values. These threshold values are renumbered to obtain the threshold values of the final 
dither threshold array. By changing the disk radius, the screen dot size can be adapted to the characteristics of par- 
ticular printing devices. Larger screen dots may improve the tone reproduction of printers having important dot 
gain. 

Keywords: Color printing, halftoning, stochastic screening, clustered-dot dithering, space-filling curves 

1. Introduction 

Traditional periodic screening based on different angles and frequencies for the different color layers is by far the 
most frequently used technique in offset color printing. 

Stochastic screening represents an alternative solution for color reproduction (see [Widmer92], [Seybold93], 
[Ostromoukhov93], [Eschbach94] , [Kang96]). It appeared in the mid 1980s and became popular in the early 1990s 
since it offers finer grain compared to traditional screening. It also permits to superimpose more than three layers 
without Moires. This method has shown its usefulness for high resolution devices at relatively high lineatures. 
However, it requires the printing process to be controlled in order to keep the dot gain within a given range of tol- 
erance. 

Today, with the advent of high-resolution desktop color printers (electro-photography, ink-jet) there is a clear need 
to re-design existing stochastic screening techniques in order to adapt them to these new devices. This is especially 
the case for electrophotographic color printers which show a non-linear behavior different from the behavior of off- 
set printers. Tuning the size of the printed dot according to the printer's characteristics may considerably improve 
the tone reproduction curve. 

The present article explains the process of building variable-size stochastic screens. We propose methods to build a 
large size dither matrix which may comprise between 2500 and 10000 stochastically laid out screen dots. In sec- 
tion 2, we describe several methods to lay out the dots within the large dither matrix. Once the screen dots are laid 
out, i.e. their centers are positioned, we show in section 3 how to create dither threshold levels ensuring a linear dot 
growth behavior across all gray levels. In section 4, we draw the conclusions. 

2. Screen dot layout within a large dither matrix 

Regarding the terminology, we use the term screen dot or simply dot for a set of black pixels clustered around a 
center point. We use the term dither cell or simply cell for a dither matrix element containing one threshold level 
and being used to generate one pixel of the output image. 



1. http://lspwww.epfl.ch/~victor 

2. Email: RD.Hersch@epfl.ch 

In Color Imaging: Device- Independent Color, Color Hardcopy, and Graphic Arts IV, SPIE Vol. 3648, pp. 496-505, 1999. Reprinted in Journal of 
Electronic Imaging, Special issue on Color Imaging: Device-Independent Color, Color Hardcopy, and Graphic Arts, October 1999. (to appear) 



The layout of the dots within a large dither matrix should respect the following requirements: 

a) Low frequencies should be avoided, i.e. the most representative frequencies present in the rendered halftone pat- 
terns should correspond to the screen dot periods. 

b) The frequencies present in the rendered halftone patterns should be isotropic, i.e. there should be no stronger 
screen frequencies at specific angles. 

c) The produced final dither matrix will cover only a part of the final image. It should therefore, when repeated hor- 
izontally and vertically, cover the whole image without creating any discontinuities. The final large dither matrix 
should therefore be continuous across its borders, i.e. the matrix should wrap around, horizontally and vertically. 
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Fig. 1. (a) Matrix covered by partially overlapping disks (dotted circles). The disks shown have a radius of r/2. Thick solid 
lines show the corresponding Delauney triangulation. Thin solid lines show the limits of the screen dots, (b) Histogram of 
screen center distances for screen dots generated with a disk radius of 16. (c) Histogram giving the distribution of screen dot 
surfaces for screen dots generated with a disk radius of 16. 



To layout the centers of the dots, we can define the minimal distance between two dots as the radius r of a disk. To 
obtain a spatial distribution of the positions of the centers of all screen dots within the dither matrix, we propose to 
first mark all cells of the dither matrix as free cells. We then define a visiting path enabling us to visit all cells of the 
dither matrix. At each visited free cell, we place a discrete disk having a radius of r cells. All cells covered by such 
a disk are marked as occupied. Since the visiting path traverses all matrix cells, the whole matrix becomes covered 
by partially overlapping disks (dotted circles in Fig. 1). 



Fourier amplitude spectrum 

dot centers, disk radius r = 8 dot centers, disk radius r = 16 of dot center images, r=16 




Fig. 2. Dot centers produced by visiting the dither array (a), (b) scanline by scanline, (d), (e) according to the Hilbert space- 
filling curve, (g), (h) according to the rotated Hilbert space-filling curve, (j), (k) according to a random space-filling curve and 
(c), (f), (i), (1) the corresponding Fourier Amplitude spectra. 



All disk centers become the center of a screen dot. To define the maximal surface of each screen dot, a Delauney 
triangulation is applied to the set of centers [Goodman97]. The corresponding Voronoi diagram may have been 
used to define the maximal surface coverage of each stochastically laid out screen dot. However, to obtain well 
clustered black dots at high intensity levels and well clustered white dots at low intensity levels, an equivalent, but 
slightly different scheme is used to create the maximal surface coverage of each screen dot (see section 3). 

2. 1 Visiting path: scanline by scanline 

The properties of the dot distribution will be determined by the visiting path. A naive visiting path scanning the 
dither matrix cell by cell and scanline by scanline produces a set of periodically laid out and tightly packed screen 
dot centers (Fig. 2a and 2b). The resulting screen is a periodic diagonally oriented screen. The basic screen fre- 
quencies are located at particular locations and are therefore not isotropic. 

2.2 Visiting path: Hilbert space-filling curve 

To improve the dot distribution, we can use as a visiting path the Hilbert space-filling curve [Sagan94]. In that case, 
the dither matrix is of size N, where N is a power of two. Since parts of the Hilbert space-filling curve incorporate 
some symmetry properties, the resulting dot center distribution is again isotropic (Fig. 2d-f). 

2.3 Visiting path: rotated Hilbert space-filling curve 

An improved distribution of dot centers is obtained by creating a visiting path made of a Hilbert space-filling curve 
rotated by a one-to-one discrete rotation [Ostromoukhov94]. The rotated space-filling path is obtained by rotating a 
multiple of the original Hilbert space-filling path with a Pythagorean angle atcin(3/4). 

Such a discrete one-to-one rotation is obtained by rotating with a Pythagorean angle a = atan(3/4) = 36.87° and by 
applying rounding operations. Let us assume that (u,v) are the coordinates of the original dither array cells (the 
lower-left cell lies at the origin) and (x,y) the coordinates of the rotated dither array cells. Then the discrete one-to- 
one rotation is expressed by: 



where c = 5, a = 4 and b = 3. 

To be horizontally and vertically repetitive, any structure such as a matrix or a discrete path to be rotated by a 
Pythagorean angle a = atan(3/4) must have a horizontal and vertical period which is a multiple of 5 
[Ostromoukhov94]. Therefore, one must choose to rotate a dither array made of a single Hilbert curve, of size 5N 
by 5N, where N by N is the size of the array containing the Hilbert curve. Fig 3 shows the original 4x4 Hilbert 
space-filling path and a rotated space-filling path obtained by rotating a multiple of the original Hilbert space-fill- 
ing path by the Pythagorean angle a = atan(3/4) = 36.87°. 

The rotated Hilbert space-filling curve is asymmetric. One can observe both in the image space (Fig. 2g and 2h) 
and in the Fourier amplitude spectrum (Fig. 2i) that the resulting dot center distribution is almost isotropic. 

2.4 Visiting path: Random space-filling curve 

Let us explore the results obtained by a random space-filling curve. A random space-filling curve where proximity 
of successive cells is not required can be constructed by starting from a list of cells, initially ordered cell by cell 
and scanline by scanline. Two cells within this ordered list are randomly selected and permuted. After a number of 
permutations, the cell visiting order is completely randomized. There is absolutely no symmetry in the sequence of 
visited cells. As expected, one can observe both in the image space (Fig. 2j and 2k) and in the Fourier amplitude 
spectrum (Fig. 21) that the resulting dot center distribution is isotropic. Furthermore, multiples of the basic fre- 
quency components are much less pronounced than in the case of the rotated Hilbert space-filling curve visiting 
path. 

Methods for finding screen dot layout presented in this article differ from other methods known in the art (see for 
example [Eschbach90] , [Marcu98]). 





Fig 3. (a) Original 4x4 Ililberl space-filling path and (b) a rolaled space-filling path obtained by rotating a multi- 
ple of the original Hilbert space-filling path by the Pythagorean angle atan(3/4). 

3. Generating the dither matrix with stochastic screen dots 

Starting from a set of stochastically distributed screen dot centers, the final dither matrix is obtained with the fol- 
lowing steps: 

a) Generate a Delauney triangulation of the set of screen dot centers 

b) From the Delauney triangulation, define the full surface coverage of each screen dot (maximal surface coverage) 

c) Generate iso-intensity lines around each screen dot for a set of different intensity levels (for example for 32 



intensity levels) 

d) With a standard software package, such as Photoshop, according to the defined iso-intensity lines, create thresh- 
old levels for all cells of the dither matrix (this corresponds to the generation of a grayscale image according to the 
given iso-intensity lines) 

e) Renumber the previously generated threshold levels so as to ensure a uniform distribution of threshold levels, 
i.e. a flat histogram of the dither threshold levels present in the dither matrix. 

Many books describe the Delauney triangulation of a planar set of dots (see [Goodman97]). Once the triangulation 
has been carried out, the full surface of each screen dot is determined geometrically by connecting triangle centers 
with their respective edge middle points (Fig. 4a). 

Iso-intensity lines are determined so as to ensure that the surface within two successive iso-intensity lines grows 
linearly with intensity. Therefore, the positions of successive intersections of iso-intensity lines and triangle bound- 
aries needs to be appropriately computed, taking into account that the quotient of the surfaces of two scaled trian- 
gles is the square of the quotient of their respective side lengths. Figure 4c illustrates the generated iso-intensity 
lines defining equally large surfaces for a certain numbers of screen dots. 




Fig. 4. Creation of a stochastic dither matrix using a rotated llilberl visiting path: (a) maximal screen dot surface 
(hatched); (b) creation of iso-intensity lines by subdividing the triangles' edges into appropriate intervals; (c) mod- 
ified iso-intensity lines; and (d) detail of the resulting dither matrix, plotted with a few discrete intensity levels 
bounded by iso-intensity curves. 



Because of the subsequent renumbering of cells, the iso-intensity lines are allowed to incorporate some variations 
and the generated surfaces are allowed to differ slightly one from another. Fig. 4d shows the iso-intensity lines gen- 
erated so as to produce screen dots with additional vertices, i.e. smoother screen dots than the screen dots produced 
in figure 4c. 

The generated iso-intensity line artwork is rasterized into an anti-aliased grayscale image using a standard imaging 
software package (for example Adobe PhotoShop). The grayscale levels are an approximation of the desired 
threshold levels and represent a surface interpolation of the previously generated iso-intensity lines. 
These grayscale levels are renumbered, i.e. identical gray levels are assigned successive numbers. The resulting 
renumbered array is the desired dither threshold array: the renumbering process ensures that dither threshold levels 
are uniformly distributed. 

The generated large size dither array is used to create halftones with a stochastic screen dot. The disk radius r 
defines the minimal distance between the screen dots. The histogram of screen dot center distances shows that 
screen dots are apart by a maximal distance of 2r (Fig. lb). Histogram lc gives a histogram of the surface associ- 
ated to screen dots, for a disk radius of 16 pixels. It shows that the minimal screen dot surface is close to that of a 
disk with radius r/2. 

Figure 4b shows the dither matrix produced by the random space-filling visiting path. To provide a sharper view of 
the screen dots, the dither matrix is plotted with a few discrete intensity levels bounded by iso-intensity curves. 
Fig 5a shows a wedge rendered with the dither matrix shown in Fig. 4b (disk radius r=16), and Fig 5b shows a 
wedge rendered with a dither matrix generated using a disk of radius r=8. Observe the generated stochastic screen 
dots: there is no predominant orientation and no low frequency patterns are visible. 




Fig 5. (a) Grayscale wedge rendered with the dither matrix shown in Fig. 4b, using the rotated I Iilbert space-lilling curve, disk 
radius r=16, printed at 600 dpi. (b) Grayscale wedge rendered with a dither matrix using the rotated I Iilbert space-lilling curve, 
disk radius r=8, printed at 600 dpi. 

Please note the non-symmetric behavior of the screen dots in highlight and shadow regions. According to the lay- 
out of the screen dots (Fig. 4b), the number of white dots is approximatively twice the number of black dots. There- 
fore, in highlight regions, black dots have a surface which is twice as large as white dots in corresponding shadow 
regions. 



4. Conclusions 

A new technique for building stochastic clustered-dot screens has been described. A large dither matrix comprising 
thousands of stochastically laid out screen dots is constructed by first laying out the screen dot centers. Screen dot 
centers are obtained by placing discrete disks of the chosen radius r at free cell locations when traversing the dither 
array cells according to either a discretely rotated Hilbert space-filling curve or a random space-filling curve. After 
Delauney triangulation of the screen dot centers, the maximal surface of each screen dot is computed and iso-inten- 
sity regions are created. This iso-intensity map is converted into an anti-aliased grayscale image, i.e. into an array 
of preliminary threshold values. These threshold values are renumbered to obtain the threshold values of the final 
dither threshold array. 

The resulting stochastic clustered dot screen incorporates screen dots whose distances vary between the chosen 
disk radius r and 2r. For the random space-filling curve, the histogram of the screen dot surface is similar to a 



Gaussian function. The largest screen dot surface is roughly 3 times larger than the smallest screen dot surface. By 
changing the disk radius, the screen dot size can be adapted to the characteristics of particular printing devices. 
Larger screen dots may improve the tone reproduction of printers having important dot gain. 
The stochastic screens produced using the presented method look like fine mezzotint (see Figs 6 and 7). The pro- 
posed stochastic screen generation method is of interest for high-resolution office printers: at 1200 dpi, a stochastic 
screen generated with a disk of radius r=8 corresponds to a traditional screen of 150 lpi. The extension of the pre- 
sented method to color images is underway. This extension will also incorporate sophisticated dot gain compensa- 
tion methods. 
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Fig. 6. A sample image, generated with a dither matrix incorporating stochastic screen dots produced with a disk of radius 16 
(random space-filling curve), printed at 600 dpi. 




Fig. 7. A sample image, generated with a dither matrix incorporating stochastic screen dots produced with a disk of radius 8 
(random space-filling curve), printed at 600 dpi. 



